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We attempt to explain the A/p ratios measured in heavy ion collisions at 11.6 A • GeV/c beam 
momentum within a hadronic framework. This ratio is enhanced relative to corresponding ratios 
in pp collisions, and is large when compared with thermal fits to heavy ion data. Using a detailed 
cascade calculation, we show that different annihilation cross-sections of A's and p's, and the net 
conversion of p's to A's, do not account for the enhancement in central collisions. For larger impact 
parameters, however, hadronic mechanisms may well suffice to produce the observed enhancement. 
Uncertainties in elementary cross-sections and formation times are considered. 



I. INTRODUCTION 



Anti-lambda (A) cross-sections in heavy ion collisions are of interest because strangeness production is a potential 
signal for the formation of a quark-gluon plasma (QGP) in heavy ion collisions jl]-f|. The ratio of anti-lambdas to 
anti-protons (p) is of particular interest since it reflects the production of strange anti-quarks relative to non-strange 
light anti-quarks, and should increase substantially relative to the production expected from a superposition of NN 
collisions, if a QGP is formed. 

Recent experiments have reported measurements of p and A production in various heavy ion systems at the 
\Q [ Brookhaven AGS and at the CERN SPS Experiment E859 reports the ratio of A and p rapidity distribu- 

tions to be 3 ± 1 ± 1 in central Si+Pb collisions Q. This ratio is corrected for finite experimental acceptance, 
efficiencies, and for p creation from A decay (feed-down), and also takes into account that neutral £ particles cannot 
be distinguished from the A sample. Experiment E864 has measured the p production cross-section in Au+Pb 
at 11.6 A • GeV/c. They compare this measurement with a similar one from the E878 collaboration that was obtained 
with a focusing spectrometer, and interpret the difference between the two measurements as an indicator of A pro- 
duction. E864 estimates that A/p > 2.3 at the 98% CL, with a most probable value of 3.5. At the SPS, NA35 has 
published A/p ratios for pp, pA, S+S, S+Ag, and S+Au collisions at 200 A • GeV/c, and they observe a significant 
rise from 0.25 for pp-collisions to 1.5 for the heavy ion systems 

It is tempting to interpret the large reported ratios as evidence for the formation of a QGP: Many authors ]9}-|ll]] 
have argued that the CERN multi-strange baryon ratios ]l2|-|l4| can only be described by a QGP scenario. Moreover, 
in a systematic study of strangeness production in heavy ion collisions from SIS to SPS energy ranges, mostly based 
on kaon ratios, Geiss et al. |15| conclude that a QGP scenario is required to explain that data at AGS energies. 
This conclusion is, however, challenged by Au+Au RQMD simulations that find that the p production cross-section 
appears to be lower than expected from the scaled NN p cross-section. 

Here we shall try to address quantitatively all possible hadronic contributions to the A/p ratio using a cascade 
model |l^] . We restrict calculations to AGS energies and show that "differential annihilation" of the two species and 
p-to-A conversion processes can indeed enhance the A/p ratio. We conclude that the effect is not large enough at 
central impact parameters, but may well suffice to explain the data from more peripheral events - thus hinting at a 
production mechanism outside of standard hadronic interactions. 

We present our arguments as follows: Section II outlines a thermal model calculation with chemical equilibration 
and feed-down contributions. We fit the parameters of the model to rather liberal ranges of various kaon, pion and 
proton ratios measured in central Au+Au events, and show that the resulting A/p ratio lies well below the E864 
value. This conclusion has also been reached by other authors |L7[[l{|. We consider transport simulations in the 
remainder of the paper, briefly describing the mode l ||l6| and discussing the relevant cross-sections in more detail. 
Actual measurements of the A and p cross-sections |20| , |2l| | are used, rather than event generator parameterizations. 
Section IV presents our results, discussing also the uncertainties in the calculations, particularly those arising from 
the relatively poorly known A annihilation in the relevant energy range, and from particle formation times. We shall 
consider first central events, and then the more recent E864 results [p2| at larger impact parameters. Speculative 
conclusions based on the quantitative discrepancy between calculations and the actual measurements at small impact 
parameters are presented in the last Section. 
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II. THERMAL MODEL 



We consider in this Section a thermal model fit for Au+Au collisions at 11.6 GeV/c, at central impact parameters 
only. Our main motivation is to show that in such events the measured A/p ratio is larger than predicted by a 
standard thermal model, even if we push the input fit ratios to extreme values. 

The model is based on thermal and chemical equilibration up to some freeze-out, so that all relative abundances 
can be obtained from four parameters: the freeze-out temperature To, the electric chemical potential fj, e , the baryonic 
chemical potential /Xf,, and the strange chemical potential [i s . Here, 

Hi = qifj, e + bifi b + Si fi s 

are particle chemical potentials, where qi, hi and Si are the charge, baryon number and strangeness of species i, 
respectively. Note that the strange chemical potential is not fixed a priori. The distribution functions we use are 
for ideal Fermions and Bosons - we do not include mean fields or excluded volume modifications p3| ], which would 
alter our conclusions quantitatively, but not qualitatively. Of course, the final abundance of a particle is determined 
by both the primary number of this species at freeze-out, and by feed-down from heavier species after freeze-out. 
We include all mesons with rest mass < 1 GeV/c 2 , and all baryons with rest mass < 1.7 GeV/c 2 . We note that our 
final A numbers throughout include S°'s, as none of the experimental measurements can distinguish between the two 
species. 

The results of the calculation are summarized in Table 1. The thermal fit parameters To, /z e , fib and fj, s are obtained 
by applying the conditions \Q/B — 0.40| < 1% and IS"! < 2%, where Q, B and S are the overall charge, baryon number 
and strangeness of the system, together with the constraints listed in Table 1. As motivated above, we purposefully 
choose rather large ranges in these latter constraints to show how difficult it would be to obtain the experimental A/p 
ratio. The "errors" on the fitted values for T and the chemical potentials listed in the caption of Table 1 indicate 
roughly the ranges that lead to results consistent with the constraint intervals .[] 

While the agreement of most ratios with the data is relatively good (see Table 1, also Refs. |l7],[l8|]), we see that 

A/p~ 1.9, at best, well below the 98% CL lower value of 2.3 reported by E864. In Figure 1 these results are illustrated 
graphically. Panel (a) shows A/p as a function of the K + /K~ ratio for various freeze-out temperatures To. Large 
A/p > 2 result only if the freeze-out temperature and/or observed K + / K~ ratio are pushed unreasonably high. A 
similar conclusion follows from Figure 1(b), which shows A/p as a function of K + /tt + . In both panels the error 
bars/scatter points indicate values consistent with constraints not shown, and the dashed line represents the 98% 
lower CL of the E864 experiment. 

Clearly, the experimental A/p ratio has, at the very least, a significant non-thermal component: the lower experi- 
mental bound reported by E864 is larger than any reasonable thermal fit would allow. 

III. CASCADE MODEL AND CROSS-SECTIONS 

Generally, a thermal description is useful for a given particle species if its mean free path is small compared to the 
system size. This fact alone means that other approaches should also be investigated. We do so in the next Section, 
and discuss here the relevant cross-sections. 



Our cascade model will be described in detail elsewhere |16|, and we give here only a brief overview. It is used 
in this work to describe the evolution of the pions, kaons, and nucleons which act as a "background" for the p and 
A species that interest us. The model is based on NN input data, and includes elastic and inelastic processes. For 
clastic collisions, resonance formation is considered if the energy and quantum numbers of the combined system are 
suitable. For inelastic collisions, we consider single diffractive processes, which effectively correspond to resonance 
production at low energies, and non-single diffractive processes, responsible for the remainder of the inelastic cross- 
section. All cross-sections are taken from experimental data whenever available, otherwise we use NN, irN or KN, 
and 7T7T or ttK cross-sections for baryon-baryon, meson-baryon, and meson-meson interactions, respectively. Average 
particle numbers, the multiplicity distribution, and species abundances are also taken from data whenever available. 
For leading particles, the fractional longitudinal momentum is determined from a uniform distribution, while the 
transverse momentum is taken from p t exp(— (3p t ) with (3 w 4.44 (GeV/c) -1 . The transverse momenta of secondary 



1 We note that these parameters are not inconsistent with overall energy conservation considerations. 
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particles are determined in a similar way, while their rapidities are chosen from a fiat distribution. All conservation 
laws are strictly enforced. 

The following processes are of particular importance for our discussion. Firstly, the production of p's and A's: 

(1) 
(2) 
(3) 
(4) 
(5) 
(6) 

Here, b represents a non-strange baryon, N a nucleon, B any baryon, M a light unflavored meson, and S ± a ±1- 
strangeness meson. We return to the question of production below. Secondly, in the same notation: 

b + p-^X (7) 

b + A^S + +n (8) 

6 + 5+^A + vr (9) 

p + M -S-A + 7T + S~ (10) 

A + M — > 6 + S + (+tt) (11) 

We refer to (0) and (||) as annihilation processes, and (^-([Tll) as conversion processes. The inverse for all 2 — ► 2 body 
processes listed is also included, but we have not considered 3, 4, ... — > 2 detailed balance. 
The experimental pp annihilation cross-section is well known, and can be parameterized as 

^ nih ( P m) = Q7p- a ° b 7 mb , (12) 

where pi a b is the momentum of the "beam" particle in GeV/c with the "target" at rest. The solid line in Figure 2 shows 
this parameterization of the data (crosses). The A annihilation cross-section, on the other hand, is relatively poorly 
known, especially in the energy range we are interested in. We model it by assuming that the elastic cross-sections 
for pA and pA are equal, and then use the data of Ref. to obtain: 

T a p r ih (Piab) = 15 ( P iab/l0)- a mb , (13) 

where the same comments apply as for Eq. (|l2|). While the data is best fit by a = 0.5, the uncertainty is rather large 
(see Figure 2; the diamonds are data from Ref. |plf ). In fact, one might well argue that the A data is consistent with 
a = 0, as shown by the dot-dashed line in Figure 2. We take a = as an extreme value that should lead to the 
greatest differential annihilation. For a — 0.7, on the other hand, the A data is practically indistinguishable from the 
p data. We shall subsequently investigate the behavior of A/p with a. 

The p and A scattering with mesons is also an important process we need to consider. Broadly speaking, we have 
three types of collision: (1) thermalization of p's and A's through elastic collisions; (2) production of resonances that 
eventually decay back into p's or A's; and (3), most importantly, net conversion of p's to A's. Chief amongst these is 

p + K + — > 7r + A (or resonances of A) , (14) 

for which we know the charge conjugate reaction to have a sizeable cross-section. The process ( |T4| ) thus contributes 
to reducing the p abundance while enhancing the A abundance in the final state. Given the pronounced strangeness 
enhancement in large systems such as Au+Au, the process ( |l4|) should be a relatively important piece of the A/p 
"puzzle." 

The production rate of p's or A's is only O(10~ 2 ) per event at the AGS, and cascade calculations are rather CPU 
intensive. We shall therefore perform an effective calculation of the survival probability by putting p's or A's in by 
hand, wherever and whenever a collision occurs in which the energy is sufficient for a pp or AA pair to be produced. 
Thus we assume that the pair production does not depend on energy, once above threshold, and take the relative 
initial production of p's and A's from pp data (see below). The presence of initial p's or A's is not permitted to 
influence the evolution of nucleons, pions and kaons - we restore particles that interacted with p's or A's to their 
pre-collision kinematics. 

We may then calculate the A/p ratio via the survival probabilities P s and the number of converted anti-lambdas 
that survive, A c : 
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- ~ r p —— + 2r c — , (15) 
p P s (p) p s 

where p s is the number of surviving p's. Here, r c rs 0.25 is a correction factor, as we overestimate A-production by 
assuming that all collisions of non-strange anti-baryons with positive strange mesons result in a A. The factor two 
accounts for the conversion of n's to A's. Finally, r p is the A/p ratio in pp collisions. At y/s ~ 20 GeV, it has a value 
p8| [n| of 0.25-0.30. At AGS energies (y/s ~ 5 GeV) its value is less established. Using Rcfs. p^-[30|], we infer a value 
of ~ 0.2, but we shall use a value of r p ss 0.25 + 0.0 —0.1 throughout this work. We note that mean held effects may 
very well reduce the effective r p for in-medium production when compared p4| to the free space values. Our chosen 
value of r p and error bars imply that the heavy ion A/p ratios we obtain in subsequent calculations will be upper 
bounds. 



IV. CASCADE CALCULATION: RESULTS 



We begin by discussing the most central events, roughly 10% of the total cross-section. To illustrate the effect of 
differential annihilation versus net conversion, we show in Table 2 the different components of Eq. ( |l5|) for different 
a and r = (no formation time) in central Au+Au simulation at AGS energies. For small a (small A annihilation 
cross-section), A/p enhancement results almost entirely from differential annihilation. As a increases, the survival 
probability of a A and a p become equal, and the sole enhancement in the hnal ratio results from conversion. Of 
course, as a increases the ratio A c /p s decreases, since A c refers to surviving anti-lambdas. 

We show in Figure 3 the hnal A/p ratios from cascade simulations (circles) for several values of a and r, and for 
various systems, all at central impact parameters. The dashed lines show the E864 98% CL values for Au+Au. We 
clearly see that reducing the difference between A and p annihilation, i.e., increasing a, reduces the ratio (Figure 3 
(a)-(c)). Similarly, a larger formation time r also leads to a smaller ratio - the effective system size for annihilation 
and conversion is reduced as all particles require some time to be formed (see panel (d)). Moreover, for Au+Au, the 
maximum ratio obtainable (for a = and no formation time) lies at about the the E864 98% CL lower limit. 

In discussing the results of Figure 3, we note that (1) realistic values for a are probably closer to 0.5 than 0.0; (2) 
formation times are probably not zero; and (3) our initial, i.e. pp A/p ratio is probably an over-estimate, particularly 
in light of possible medium effects that penalize A production relative to p production. This latter fact is reflected in 
the error bars for the ratio - our choice for r p = 0.25 +0.0 —0.1 is what we consider an upper bound. We conclude 

that for central Au+Au events the most likely value of A/p is most definitely ~ 2, and more likely somewhat below 
unity - smaller than the 98% CL reported by E864. 

The A/p ratio is larger in Si+Pb collisions than in Au+Au collisions in our calculation, but it is not a large effect 
(compare Figure 3 (a) and (c)). Within our model this difference has a simple geometric explanation: For Si+Pb, the 
p's and A's are produced close to the beam axis, leading to a larger average nuclear thickness that the anti-particle 
must traverse. Once again, though, realistic input parameter values of r ~ 1 and a ~ 0.4 seem to suggest that our 
hadronic cascade misses the experimental lower bound of E859. 

The solid lines in Figure 3 are results from a model calculation [p2[ that uses annihilation and net conversion mean 
free paths as an input to rate equations in a simplified geometric picture. Its only free parameter mocks up the effect 
of local matter expansion, local average momentum distributions, etc. It reproduces the trends in the data rather 
well. 

Recently, E864 has also reported the impact parameter dependence of the A/p ratio in Au+Au |22]| . In Figure 4, 
the upper and lower solid line (drawn to guide the eyes) show the most probable (diamonds) and 98% CL (squares) 
values of these results, respectively. There are four centrality bins: From left to right they are < 10%, 10-30%, 
30-70%, and > 70% of the total cross-section, as obtained from transverse energy cuts. The short-dashed line joins 
cascade calculations (crosses) for a = and r = 0. Recall that both these parameters values are extreme: a = 0.4 
(box-crosses, joined by long dashes in Figure 4) is more likely, and does not reproduce the central data bin result. For 
larger impact parameters, however, we see that this conclusion is considerably weakened - the cascade result initially 
rises with increasing impact parameter to far above the data for a = 0, while for a = 0.4 it remains relatively flat, 
but eventually does intercept the experimental data in semi-central collisions. 

Thus, while hadronic mechanisms seem inadequate for the most central events, we are not able to conclude with 
certainty that they fail to describe the data at larger impact parameters. We note that at higher energies the difference 
between the A and p annihilation cross-section practically vanishes. Within the differential annihilation model, one 
might therefore predict (in hindsight) that A/p at the SPS should be lower than at the AGS, as indeed observed, and 
indicating that the model has the correct beam-energy dependence. 
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V. CONCLUSIONS 



The A/p ratio in NN collisions at AGS energies is ~0.25, heavy ion experiments measure much larger values. 
A hadronic thermal model with chemical equilibration and feed-down contributions can only be pushed to values 
above 2 with much difficulty For Au+Au, using the (unreasonable) extremes of K + /n + ~ 0.30, K + / K~ <~ 6.0, and 
T <~ 140 MeV, we obtain an upper limit of A/p ~ 2, below the E864 98%-lower confidence limit of 2.3. 

We have considered a non-equilibrium description, and presented results from a detailed cascade calculation. Chief 
inputs to our calculations are the A and p annihilation cross-sections with nucleons, as well as processes that convert 
p's into A's. The A annihilation cross-section is relatively poorly known in the relevant energy range, but it may 
be argued that it is somewhat less than p annihilation at AGS energies. Our parameter a controls this difference, 
with a = corresponding to the greatest difference, and a = 0.7 corresponding to a cross-section that is practically 
indistinguishable from p annihilation. Currently, the cross-section data is ("best") fit by a — 0.5. Another important 
input is the cross section of NN A to p production, which in our calculation may well be an overestimate, compared 
to free AW data and also in view of medium effects on the particle production. Our asymmetric error bars on the 
input AW ratio reflect this recognition. 

In Au+Au cascade simulations for small impact parameters at AGS energies, the largest A/p ratio obtainable is 
<~ 2.4+0.2 — 1.0, for a flat and small A annihilation cross-section (a = 0) and no formation time (r = 0). Both these 

input values are extreme. What we currently believe to be more reasonable inputs lead to a value of A/p~l (~ I 
for Si+Pb), far below the lower bound of the E864 experiment. The model essentially fails to describe the ratio. For 
larger impact parameters, however, we cannot conclude with certainty that the hadronic mechanisms are inadequate 
in describing the data. While the data drops dramatically as the impact parameter increases, the calculation remains 
relatively flatter, and even rises somewhat for an energy independent A annihilation cross-section. 

In conclusion, the large A/p ratios in central Au+Au collisions at 11.6 A • GeV/c are not easily explained by hadronic 
mechanisms. This is no longer true at larger impact parameters, and it is interesting to speculate if QGP formation 
in central events might be a possible explanation. Of course, more accurate measurements of the A annihilation 
cross-section in the relevant energy range are needed, as well as a better understanding of hadronic medium effects 
in A and p production and annihilation, before any more definite conclusion can be reached. 
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VII. TABLES 



TABLE 1. Thermal and experimental particle ratios for central Au+Au collisions at 11.6 GeV/c. The parameter 
ranges are T = 120 ±14 MeV, ^ b = 556 ±19 MeV, ^ s = 111 ± 14 MeV, and [i e = -14±2 MeV. The second column 
lists the range of the various constraint ratios considered. 



Ratio 



A-+/7T+ 

K+/K- 

K-/TT- 
7T + /p 

n~/p 
A/p 
p/p 
A/p 



Constraint 
Range 



Thermal Model 



0.23 ±0.03 
4.73 ±0.53 

(3.50 ±0.62) x 10~ 2 
0.71 ±0.09 
1.00 ±0.10 
0.16 ±0.02 

(3.48 ±3.44) x 10~ 4 
1.58 ±0.30 



Data 



Ratio 



Rapidity 



Ref. 

~|2f| 

[25] 



0.16-0.28 
4.0-6.0 

0.6-1.2 
0.8-1.4 



0.22 ±0.01 
5.0 ± 1.0 
0.028 

1.00 



0.5-1.3 
0.5-1.3 
1.2-2.0 

1.2-2.0 



(5 ±2) x 10~ 4 
3.5 



1.6 (p±~0) 
1.6 (p±~0) 
> 2.3 (98% CL) 1.6 (p_L~0) 



TABLE 2 Survival probabilities P s , the ratio of surviving converted A c 's to surviving p s , and final A/p ratio for 
various a and r = fm/c, in central (~10%) Au+Au cascade simulations. 



a 


P S (A) (%) P,(p) (%) A c /p s 


A/p 


0.0 
0.2 
0.4 
0.7 


12.1 ±0.3 1.5 ±0.1 0.75 ±0.07 
6.0 ±0.2 1.5 ±0.1 0.45 ±0.06 
3.3 ±0.2 1.5 ±0.1 0.38 ±0.06 
1.5 ±0.05 1.5 ±0.1 0.16 ±0.02 


2.4+0.2-1.0 
1.2 + 0.1-0.5 
0.7+0.1-0.3 
0.3+0.04-0.1 
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VIII. FIGURE CAPTIONS 



FIGURE 1. (a) Thermal A/p ratios as a function of K + /K~ for various freeze-out temperatures To. (b) Thermal 
A/p ratios as a function of K + /tt + . The error bars/scatter points indicate values consistent with constraints not 
shown (see Table 1 for details). The solid lines are to guide the eye only; the dashed line represents the 98% lower 
CL of the E864 measurements. 



FIGURE 2. The annihilation cross-sections of p (A) with nucleons, as a function of the incident momentum of the p 
(A). The solid and dashed lines are fits to the p and A data, respectively, while the dot-dashed curve is the extreme 
fit to the A data (with a = in Eq. (|l3|)). The A data (diamonds) are from Ref. pl| . 



FIGURE 3. A/p ratios calculated in the cascade model, for various values of a, and t for central Si+Pb and Au+Au 
collisions. The solid lines are results of a geometric model calculation of Ref. |32|, while the dashed lines show the 
E864 lower 98% CL 101. 



FIGURE 4. A/p ratios as a function of centrality for Au+Au and zero formation time r. All lines are drawn to 
guide the eye only. The upper solid line shows the E864 most probable values, the lower solid line represents the 
E864 98% CL ||^2|]. The crosses are cascade calculations for a — 0, while the cross-box points are for a = 0.4. Some 
points are offset slightly to clarify the error bars. 
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